Spontaneous breaking of the Fermi surface symmetry in the t — J modeh a numerical 

study 
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We present a variational Monte Carlo (VMC) study of spontaneous Fermi surface symmetry 
breaking in the t — J model. We find that the variational energy of a Gutzwiller projected Fermi 
sea is lowered by allowing for a finite asymmetry between the x- and the j/-directions. However, 
the best variational state remains a pure superconducting state with d-wave symmetry, as long 
as the underlying lattice is isotropic. Our VMC results are in good overall agreement with slave 
boson mean field theory (SBMFT) and renormalized mean field theory (RMFT), although apparent 
discrepancies do show up in the half-filled limit, revealing some limitations of mean field theories. 
VMC and complementary RMFT calculations also confirm the SBMFT predictions that many-body 
interactions can enhance any anisotropy in the underlying crystal lattice. Thus, our results may be 
of consequence for the description of strongly correlated superconductors with an anisotropic lattice 
structure. 

PACS numbers: 74.20.Mn, 71.10.Li, 71.10.Fd 



I. INTRODUCTION 

The phenomenon of high temperature superconductiv- 
ity has led to intensive debates about possible supercon- 
ducting states in purely repulsive models, such as the 
single band Hubbard model in two dimensions. Var- 
ious strong coupling approaches show that a repulsive 
onsite (Hubbard) interaction U may lead to instabili- 
ties of the two-dimensional (2D) Fermi sea, such as an- 
tiferromagnetism, d-wave superconductivity etc. An in- 
teresting possibility is the tetragonal symmetry break- 
ing of the 2D Fermi surface due to strong electron cor- 
relations, which may result in a quasi one-dimensional 
(ID) state. This instability, which competes with, and 
seems to be suppressed by the d-wave pairing state, was 
first reported by Yamase and Kohno within slave bo- 
son mean field theory (SBMFT) in the t — J model. ^ 
The same effect was seen within a renormalization group 
(RG) study in the Hubbard model by Halboth and Met- 
zner who called it the "Pomeranchuk instability" (PI) of 
the Fermi surface. Subsequently, several authors studied 
this problem.^'^'-'^-'^i'^i^iQd°i"d^i^3d'^i^^ It was also argued 
that the tendency towards a quasi ID state may enhance 
a bare anisotropy of the underlying lattice structure, as 
for example present in YBa2Cu306+j, (YBCO).^^ 

Previous studies of the PI for the t—J and the Hubbard 
models were mainly based on SBMFT and RG methods. 
Clearly, it is desirable to use an alternative approach, 
and examine if indeed a PI is present. In this paper, 
we analyze this problem comprehensively using a varia- 
tional Monte Carlo (VMC) methodJ^-iiiSiia Our VMC 
results agree qualitatively with previous SBMFT studies 
over a wide doping range except for discrepancies that 
show up very close to half-filling. We also supplement 
the VMC results by presenting results from renormal- 
ized mean field theory (RMFT)2S4Si calculations. These 



results agree well with those from SBMFT, but reveal 
similar deviations from the VMC results near half-filling. 
We argue that the disagreement stems from limitations 
of the mean field theories when dealing with nearly half- 
filled states. 

The outline of the paper is as follows: In Section II, we 
introduce the t — J Hamiltonian and describe the VMC 
scheme. The VMC results for a quasi ID state in an 
isotropic t — J model are presented in Section III. We 
discuss the magnitude of asymmetry, the condensation 
energy per site and the model parameter dependence as 
a function of the hole concentration (doping). In Sec- 
tion IV, we provide an RMFT study for the PI and give 
an explicit comparison of the Gutzwiller renormalization 
scheme (GRS) with the VMC results. Section V is dedi- 
cated to a discussion of RMFT and VMC results for the 
anisotropic t—J model and is followed by our conclusions. 



II. MODEL AND VMC SCHEME 



We consider the t — J model in two dimensions, 



H = - ^trc]„c 



S,; Sj - -n,nj ] , (1) 



defined in a projected Hilbert space where double occu- 
pancies are forbidden. In the above equation. Si and rij 
are respectively, electron spin and density operators at 
site i. The kinetic energy is determined by t-r, the rth 
neighbor hopping integral. Throughout this paper, we 
denote nearest neighbor (n.n.) hopping terms by t, and 
choose next n.n. hopping terms t' — oi t' = — 0.3i. For 
simplicity, we only consider the superexchange coupling 
J — At'^/U between n.n. sites, and choose the value, 
J = OM. 
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We implement the no-double occupancy constraint in 
the t — J Hamiltonian by working with Gutzwiller pro- 
jected states 1 4*) = Pg |*o), where the Gutzwiller pro- 
jection operator, Pq ~ J^j(l — ni-^rii^). Such wavefunc- 
tions were initially proposed as variational states to de- 
scribe superconductivity in the proximity of a Mott insu- 
lating phase .iSiiLSSiSSi The VMC technique which allows 
for a numerical evaluation of expectation values in these 
states successfully predicted d— wave pairing in the t — J 
modeliiSiii The method has been extended more recently 
to study the coexistence of (and competition between) 
various phases like antiferromagnetism, flux states, and 
superconductivity in the t — J model^'^i^^i^^i^^. Moti- 
vated by the phenomenology of the high temperature su- 
perconductors, other improvements such as the inclusion 
of longer range hopping terms^^, increase of the num- 
ber of variational parameters, and introduction of long 
range correlations through Jastrow factors^^ have also 
been proposed. Here, we extend previous works by al- 
lowing a symmetry breaking between the x— and the 
y— direction in the variational wavefunction \^). 

The ground state wave function is written in the form, 
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FIG. 1: (color online) (a) VMC results for condensation en- 
ergies per site Ccond of the quasi ID state (A^ = 0) and the 
d-wave state (A^ 7^ 0) with t' — —0.3t. Optimal variational 
parameters Aifl and t'^^^ of the d-wave state are shown in 
(b) and (c). The errors in (b) and (c) are AAIS = 0.05 and 
Atvar = 0.05, respectively. System sizes: L = 11^ + 1 = 122, 



13^ + 1 = 170, L : 



226, and L = 17^ + 1 = 290. 



where I^'bcs) = {j2 
electron BCS wave function with 



= ^g|*bcs) , (2) 

N/2 



j |0) is the A^- 
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(a) RMFT (7t,0) 



(3t,7t) (b) VMC: L=15, x=0.07 
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To allow for a possible quasi ID state as well as for finite 
d- and s-wave pairing, we choose, 

Cfc = - 2[(H-Cr)cOSfc, + (l-C)cOsfc,] 

- 4t(,3^j. cos /Ca; COS fcy - /ivar (4) 



FIG. 2: Fermi surface of the isotropic t — J model with J = 
0.3t and t' = -0.3t at 2; = 0.07 (a) RMFT results for the 
Fermi surface of the normal state with Afc = (solid line) 
and the optimal d-wave state (dashed line), (b) Best quasi 
ID state on a (15^ -|- l)-sites lattice by VMC; filled circles 
indicate the Fermi surface. 



and 

Afc = A(,^^|,(cos - COS ky) A(,f,|,(cos fc^ + cos ky) . (5) 

We then have the following five variational parameters: 
(i) the asymmetry S^^^. between n.n. x— and y— hopping 
matrix elements; (ii) the variational next n.n. hopping 
term t'^^^; (iii) a variational chemical potential /ivar; (iv) 
and (v) variational parameters for d— and s— wave pair- 
ing, Avar and Aiar, respectively. Using standard VMC 
techniques^S, we compute energy expectation values and 
minimize the energy by searching for the optimal set of 
variational parameters. We use tilted lattices with pe- 
riodic boundary conditions. Typically, 10^ Monte Carlo 
steps are performed for each set of variational param- 
eters. The resulting statistical errors are given in the 
relevant figures by error bars. 



III. VMC RESULTS FOR THE ISOTROPIC t - J 
MODEL 

We first present results for the isotropic t~ J Hamilto- 
nian, with model parameters J = O.M and t' — — 0.3i. 
These are reasonable model parameters for the phe- 
nomenology of the high temperature superconductors. 
The optimal solution for various values of hole concen- 
tration, X = — 0.3, is determined by searching in the 
whole variational parameter space. We find that the pure 
isotropic projected d— wave state always optimizes the 

( s) 

ground state energy, i.e., the s— wave parameter Avar and 
the asymmetry vanish for all values of x, within our 

numerical resolution [AAvil = 0.05 and AJ-J^. — 0.05]. 

In Fig.^a) (circles) we show the condensation energy 
per site, econd, of the optimal state with respect to the 
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FIG. 3: (color online) VMC results for a quasi ID state for 
the isotropic t ~ J model (J = 0.3t) with t' — (circles) 
and t' = — 0.3t (squares). Doping dependence x of (a) the 
condensation energy per site, econd, and (b),(c) the optimal 
variational Sl^^. The errors in (b) and (c) are A5l° =0.05. 
System sizes: L = 15^ + 1 = 226 and L = 17^ + 1 = 290. 
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FIG. 4: (color online) VMC results for the quasi ID state in 
the isotropic t — J model (J — 0.3t) with t' — (circles) and 
t' — —0.3t (squares). Doping dependence x of the optimal 
variational t'^^,.. At^ar =0.1; System sizes: L = 15^ + 1 = 226 
and L = 17^ + 1 = 290. 



projected isotropic Fermi sea. The condensation energy, 
Gcond, is calculated by comparing the VMC energy ex- 
pectation values in the projected Fermi sea and optimal 
d— wave states. We see a continuous increase of | econd I 
and of the superconducting d— wave parameter [shown in 
Fig.^b)] as doping x decreases. The optimal variational 
value for t'^^^. is given in Fig. QJc). These results show 

that the additional variational parameters aIII- and 
are not relevant for improving the ground state in the 
isotropic t — J model. 

To uncover the PI, it is necessary to suppress super- 
conductivity by setting Afc = 0. Doing so, our VMC cal- 
culations indeed reveal a PI; we obtain an improvement 
of the energy expectation value relative to the isotropic 
projected Fermi sea by using a finite asymmetry <5^^., al- 
though the underlying lattice is still isotropic. We com- 
pare the condensation energy of this state to that of the 
d— wave state. Results are shown in Fig. Qla). As for 



the d— wave state, |econd| for the quasi ID state initially 
increases as x decreases. However its energy gain is much 
less than that of the d— wave, and so the latter is always 
favored on an isotropic lattice. Furthermore, note that 
the condensation of the quasi ID state saturates and fi- 
nally vanishes very close to half-filling. We will come 
back to this point later. 

The resulting VMC Fermi surface of the optimal quasi 
ID state at a; = 0.07 is shown in Fig. [2Ib). It reveals 
why finite size effects become important in the VMC cal- 
culations when dealing with a finite asymmetry in the 
projected Fermi sea. Varying (5^^ causes discontinuous 
changes of the Fermi surface on a finite lattice. The oc- 
cupied states regroup for certain leading to small 
yet discontinuous changes in the FS as a function of the 
variational parameter (5^^,. 

These finite size effects cause a rather large error for 
the optimal value of asymmetry, AS^^^ — 0.05, and for 
the effective next n.n. hopping, At'^^j. — 0.05 — 0.1. The 
jump size increases with decreasing system size, thus re- 
quiring sufficiently large lattices. The problem is less 
severe when considering a superconducting state, where 
the occupancy in momentum space changes continuously 
at the Fermi surface. 

To consider the effect of t' on the PI, we compared 
the two cases, t' = and t' — — 0.3t, in the absence of 
superconducting order (Afc = 0). The PI is stronger for 
t' < 0, yielding a larger condensation energy [Fig.OIa)]. 
As seen in Fig. |SIb), the anisotropy is significant even 
at higher doping levels and exists up to a; ~ 0.2 for 
t' — —0.3t. For t' = 0, the (5^^ is significant only in the 
range, x = 0.03 - 0.10 [Fig. E^c)]. In Fig. H we show the 
optimal variational value for t'^^^. Interestingly, t'^^^ 
for X ^ even for a bare dispersion t' = —0.3t. Re- 
cently, we reported a similar renormalization of the next 
n.n. hopping terms due to strong coupling effects within 
RMFTiSLffi 

Although there is good overall agreement between our 
VMC data and SBMFT results, we find clear and sig- 
nificant discrepancies in the limit x — *■ 0. As seen in 
Fig-Eta)-(c) the asymmetry goes to zero at a; = within 
our VMC calculations. On the other hand, SBMFT as 
well as RMFT (discussed in the next section) predict a 
pure ID state at half-filling when Afc = 0. This hints 
at limitations of the mean field theories when treating 
states near half-filling. We shall discuss this in more de- 
tail, after discussing results from the RMFT calculation 
for the PI. 



IV. RMFT AND GUTZWILLER 
RENORMALIZATION FOR THE ISOTROPIC 
t J MODEL 

In this section, we present some results from renormal- 
ized mean field theory (RMFT). The RMFT gap equa- 
tions have been derived in previous works^SiSi. Here, 
we solve them allowing for possible anisotropic solutions. 
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FIG. 5: (color online) RMFT results for the isotropic t — 
J model with J = 0.3t. (a) Condensation energy per site 
relative to the isotropic non-SC solution (Afe = 0, = Cm) for 
the quasi ID state (A^ = 0) and for the optimal d-wave state 
with t' = — 0.3t. (b) Asymmetry and (c) condensation energy 
per site for the quasi ID state with t' — and t' = —0.3, 
respectively. 

We find that the optimal self consistent state remains 
a pure d-wave superconductor. This is consistent with 
the VMC data presented in the previous section, and 
with SBMFT calculations^ The condensation energy of 
this optimal state relative to the isotropic renormalized 
Fermi sea is shown in Fig. |S1 To uncover the PI, we now 
look for solutions with the constraint = 0, as before. 
Doing so, we find two self consistent solutions for suf- 
ficiently small doping: an isotropic and an anisotropic 
renormalized Fermi sea. We compare the condensation 
energy of the anisotropic quasi ID state with that of the 
d-wave pairing state for a bare dispersion t' = —0.3t. 
As seen from Fig. EJa), the doping dependence as well 
as the magnitude of the condensation energy agree well 
with the VMC results in Fig. da). However, the two 
results begin to differ in the vicinity of a; = 0, where the 
condensation energy of the quasi ID state vanishes in the 
VMC scheme. 

Within RMFT, the order parameter characterizing the 
asymmetry, 5rmft is given by, 

rlD ^ - ty 
"RMFT ~ X ,X ■ 

Here, the effective hopping in the r-direction is related 
to the bare hopping tr bjtSi^ 

ir = gtU + ^^''g J £,T , (7) 

where J = 4t^/[/, = J2a(4,a(^i+r,a)\yS/o)^ and T^x,y. 
The renormalization factors for the kinetic energy, gt — 
2x/{l + x), and for the spin-spin correlation, gs = 4/(1 -f 
x)'^ , are derived within the Gutzwiller Renormalization 
Scheme (GRS). 



The origin of the PI can be seen directly from Eq. j?)). 
It shows that the tendency to a quasi ID state stems 
from the J-tcrm [2nd term in Eq. {Tj)] because it in- 
cludes the factor ^t-. For an isotropic dispersion it is the 
only quantity in Eq. ((TJ that may cause an anisotropy 
in the effective hopping tr- Similar arguments apply for 
the enhancement of a bare asymmetry 5q^ in a slightly 
anisotropic lattice (see also discussion below). The ori- 
gin of the PI may also be understood in the framework 
of a Landau-Ginzburg analysis as shown by Yamase and 
Kohno within SBMFT. ^ 

In Fig. [SJ^b) and (c), we show our results for 5^mft 
and the condensation energy relative to the projected 
(symmetric) Fermi sea. Results are shown as a function 
of X for two values of t' , t' — and t' — —0.3t. We 
see that a negative t' favors the PI, in agreement with 
SBMFT and VMC results. The asymmetry 5rmft 
shown in Fig. \^h) and agrees quantitatively with the 
VMC results in Fig. |3^b),(c) for x > 0.07. However, for 
a; — > the asymmetry increases strongly within RMFT, 
whereas it saturates and finally disappears within the 
VMC scheme. 

The Fermi surface of the quasi ID state from RMFT 
at X = 0.07 is shown in Fig. [2Ia) and agrees well with 
that obtained from VMC calculations for the same dop- 
ing [Fig.EJb)]. 

We now turn our attention to the discrepancy between 
VMC and RMFT (and SBMFT) results near half-fiUing. 
In Fig.|H| we plot the kinetic energy E^kin and the superex- 
change energy Ej (J-term in the Hamiltonian) for the 
projected Fermi Sea (A^ = 0) as a function of the asym- 
metry (5yfJ for various doping concentrations. We com- 
pare the VMC data with the results from the Gutzwiller 
renormalization scheme (GRS) for t' = 0. Fig. shows 
that the GRS only approximately agrees with the nearly 
exact VMC results. In particular, the behavior of the su- 
perexchange energy as a function of is qualitatively 
different in the two schemes. The discrepancies are acute 
near half filling (see data for x = and x = 0.035). A 
reasonable agreement is obtained for larger doping. In 
general, the VMC data shows a weaker dependence of 
the energy on the asymmetry 6^]^^. 

Our VMC results are also consistent with previous 
studies of the Gutzwiller wavefunction at half-filling in 
the ID limiljS^, which corresponds to (5^^ = 1 in our 
calculations. The superexchange energy (deduced from 
the n.n. spin-spin correlations) of the pure ID state on 
an isotropic 2D lattice is much worse than that of the 
projected isotropic 2D Fermi sea. On the contrary a 
simple evaluation by the GRS would favor a pure ID 
state at half-filling. We thus think that the drawbacks 
of the mean field theories partially stem from invoking 
ID physics by allowing a finite asymmetry. We further 
note that the renormalization effects become largest near 
half-filling. Therefore, it is naturally that discrepan- 
cies can show up in this limit, as already seen in pre- 
vious VMC studies, e.g., for the quasiparticle weight 
renormalization^^i^'^1^^1^^ . With that in mind, our VMC 
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FIG. 6: (color online) Comparison of the Gutzwiller renormal- 
ization scheme (GRS, solid lines) and VMC results (dashed 
lines with symbols) for (a) the kinetic and (b) the superex- 
change energy of the projected Fermi Sea (A^ = 0) at differ- 
ent doping levels x. The dependence on the asymmetry 5^^? 
is illustrated. VMC results are taken from the L = 15^ + 1- 
system and t' = 0. Statistical VMC errors are much smaller 
than the symbol size. 



results are not really surprising. Though the GRS gives 
valuable insights into the behavior of strongly correlated 
electronic systems, the above results show that a verifica- 
tion of GRS results by the VMC technique is often very 
important. These limitations of the GRS directly impact 
the RMFT (which is based on the GRS) and the SBMFT 
(which is closely related to RMFT). 
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FIG. 7: RMFT and VMC resuhs for the d + s-wave ground 
state of the anisotropic t — J model with J — 0.3t and 5]^ = 
{tx — ty)/{tx + ty) = 0.05. (a) Effective asymmetry 5^^ = 
(tx — ty) / {tx + ty) from RMFT as a function of hole doping 
X for (dashed) t' = and (solid) t' = -0.3t. VMC results 
for t' = are given by squares and circles for L — 122 and 
L = 170, respectively, (b) RMFT Fermi surface (solid lines) 
of the d-l- s-wave ground state and the tight binding dispersion 
(dashed) at a; = 0.08 with t' = ~0.3t and 5?d = 0.025. 



contributionr^ Interestingly, the bare anisotropy S}^ of 
the lattice is enhanced due to the electron correlations. 
Here we re-examine this prediction within the RMFT 
and VMC schemes. Results from RMFT are shown 
in Fig. [7fa) and (b) and agree quantitatively with the 
SBMFT data. As seen in Fig. |3a), the bare asymmetry 
of (5g^ = 0.05 increases within the RMFT calculations up 
to about iJ^pt ~ 0.2 in the underdoped regime. These re- 
sults are confirmed to some extent by VMC calculations 
for t' = in Fig.|7{a) (circles and squares), that show an 
increase of the asymmetry up to about (5^^ « 0.1. How- 
ever, owing to numerical difficulties, the errors in these 
VMC calculations are quite large^^^ In Fig.[7fb), we com- 
pare the Fermi surface obtained from the bare dispersion 
(i5q^ = 0.05) with that of the optimal superconducting 
state obtained by solving the RMFT equations self con- 
sistently, for a: = 0.08. As seen in the figure, the enhance- 
ment of anisotropy due to strong correlations may even 
lead to a change in the topology of the underlying Fermi 
surface. 



VI. CONCLUSIONS 



V. RMFT AND VMC CALCULATIONS FOR 
THE ANISOTROPIC t - J MODEL 

Our results from VMC and RMFT confirm that a quasi 
ID state is always suppressed by the d-wavc pairing state. 
A PI occurs only when = 0. However, the situa- 
tion can be quite different when the underlying lattice 
structure is anisotropic. In this case, the tendency to- 
wards a quasi ID state is present even in the supercon- 
ducting state. SBMFT^ predicts an optimal state which 
has a dominant d-wave symmetry with a small s-wave 



In this paper, we performed a Variational Monte 
Carlo (VMC) study of the Pomeranchuk instability in 
the isotropic t — J model. We also supplemented the 
study with results from renormalized mean field theory 
(RMFT) . Our results are in agreement with earlier calcu- 
lations from slave boson mean field theory (SBMFT) and 
exact diagonaliziation on small clusters^ and show that 
the instability is uncovered when the (optimal) d-wave 
superconducting state is suppressed. The Pomeranchuk 
instability is seen for values of hole concentration between 
X « 0.03 and x ~ 0.20, depending on the magnitude of 
the next n.n. hopping integral t' . We find that a small 
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bare asymmetry is enhanced within an anisotropic 
t — J model even in the superconducting state. Ahhough 
there is good overall agreement between the VMC calcu- 
lations and results from RMFT and SBMFT, discrepan- 
cies arise very close to half-filling, showing the limitations 
of these mean field theories in this regime. This under- 
lines the importance of complementary VMC studies to 
confirm or disprove the results from RMFT or SBMFT 
calculations. 

The tendency towards a quasi one dimensional state 
in a strongly correlated electron systems is mainly gov- 
erned by the superexchangc J. Since J oc At'^/U, a small 
asymmetry in the bare hopping integral t becomes twice 
as large in the superexchange energy. Hence, it is natu- 



ral that the effects discussed in this paper are largest in 
the underdoped regime, where the dispersion is mainly 
determined by J. The tendency towards a quasi one 
dimensional state may be also enhanced if phonons are 
coupled to the lattice. The influence of such dynamic 
effects should be clarified in future studies. 
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